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Abstract 

The influence of the material texture (substructure) on the force driv- 
ing the crack tip in complex materials admitting Ginzburg-Landau-likc 
energies is analyzed in a three-dimensional continuum setting. The the- 
ory proposed accounts for finite deformations and general coarse-grained 
order parameters. A modified expression of the J-integral is obtained to- 
gether with other path-integrals which are necessary to treat cases where 
the process zone around the tip has finite size. The results can be ap- 
plied to a wide class of material substructures. As examples, cracks in 
ferroelectrics and in materials with strain-gradient effects are discussed: 
in these cases the specializations of the general results fit reasonably ex- 
perimental data. 

1 Introduction 

To analyze the behavior of cracks, one needs to understand how the interactions 
in the material cooperate to drive the crack tip. In the common setting of 
dcformable simple bodies (Cauchy's model) the question has been clarified in 
many aspects from theoretical and computational points of view . It is not 
so for complex bodies where experiments suggest that the influence of material 
texture (substructure) on the expression of the 'force' driving the crack tip may 
be prominent. For instance, for materials that fail by decohesion or cleavage 

'From the pioneer works of Griffith (1920), Atkinson & Eshelby (1968), Baremblatt (1972), 
Preund (1972), Rice (1968), to contributions like (Adda-Bedia, Arias, Ben Amar & Lund, 1999; 
Dolbow, Moes & Belytschko, 2001; Freund, 1990; Gurtin & Shvartsman, 1997; Heino & Kaski, 
1997; Moes & Belytschko, 2002; Obrezanova, Movchan & Willis, 2002; Oleaga, 2003; Slepyan, 
2002). 



at the atomic scale, the predictions of the standard theory of fracture are not 
satisfactory 2 , and the same inaccuracy occurs in the case of cracks propagating 
along metal-ceramic interfaces. Moreover, e.g., the spontaneous polarization in 
ferroelectrics influences rather stronghly the crack growth 3 . 

The investigation is matter of continuum theories more intricated than the 
standard elasticity where only some global aspects of the material texture are 
considered through material symmetries, and a direct description of the config- 
uration of the substructure (together with its changes) and the related interac- 
tions is absent. 

However, if one would like to determine the appropriate expression of the 
force driving the crack tip in each special case of complex material, one would 
have to construct only a zoo of possible modifications of the standard theory of 
fracture for small and large strains, obscuring the basic fact that there exists a 
common unifying physical mechanism. The aim of the present paper is just to 
show such a mechanism common to all complex materials admitting Ginzburg- 
Landau-like energies and to provide accordingly a general expression of the 
driving force which can be specified and tested in special cases. We consider 
large strains to encompass cases in which they are prominent as, e.g., in presence 
of elastomers. 

In what follows, an order parameter v is assigned to each material element 
as coarse grained geometrical descriptor of the substructure. To maintain gen- 
erality, we require only that v belongs to a differcntiablc manifold M.. Each 
special choice of Ai characterizes the model of each complex material. Interac- 
tions developing extra power on the rate of the order parameter, i.e., on the rate 
of change of the substructure, are considered and satisfy appropriate balance 
equations. Models of fiber-reinforced composites, ferroelectric and magnetoelas- 
tic solids, interacting elastomers, microcracked and multiphase solids fall, e.g., 
within such a general approach. 

In what follows, first we derive balance equations at the crack tip for the 
interactions generated by the substructure, considering possible substructural 
inertia effects at the tip. Then we determine the contribution of substructural 
interactions to the expression of the driving force at the tip of the crack. It 
results in a modified expression of the J-intcgral: 

J=n 'X P (p(\\\ ± \\ 2+k ^ i/ )) i+p ) n > 

where n is the direction of propagation of the crack, the second-order tensor F is 
given by P = xpl — F T T— Vv T *S, with ip the free energy density, F the gradient 
of deformation, T the first Piola-Kirchhoff stress, <S a measure of substructural 
interactions called microstress, I the unit tensor, p the density of mass, while the 
other terms account for standard and substructural kinetic energy (if existent). 
J tip indicates a special limit process which consists in evaluating an integral on 

2 See the remarks about the interpretation of experimental results in (Wei & Hutchinson, 
1997). 

3 Experimental data and detailed analyses arc collected in (Fulton & Gao, 2001). 
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the boundary of a disc centered at the tip in a plane orthogonal to the tangent 
of the tip at a given point and in shrinking the disc up to the tip. The product 
* is defined below. Special expressions for the J-integral follow once the order 
parameter is specified. Here, to suggest examples, ferroelectrics and materials 
with strain-gradient effects are considered: the relevant J-integrals provide val- 
ues of the driving force close to experimental data. Other path-integrals besides 
J-integral are obtained: they allow us to express the energy dissipated during 
the evolution of the crack in the case in which the process zone around the tip 
has finite size. In the expression of J, two elements mark the difference with 
the standard theory of fracture, namely the densities pk (u, u) and Vu T ±S that 
underline the explicit influence of the material substructure. Whilst the former 
is negligible unless the substructure oscillates at very high frequencies, the lat- 
ter may be crucial to justify the discrepancies between experimental data and 
the previsions of the standard fracture mechanics. (To render the formulas as 
concise as possible, in the integrals we do not write explicitly line, area and 
volume differentials: the kind of integration is clear looking at the domain of 
integration directly 4 .) 



2 Order parameters for the substructure and 
the geometry of the crack 

Let B be the regular region of the three-dimensional Euclidean point space £ 3 
occupied by the body in its reference place. A generic point of B is indicated 
with X. A 'standard' deformation is described by a sufficiently smooth bijection 
89X i — ► x (X) e £ 3 (with the current place x (B) a regular region) which is 
also orientation-preserving in the sense that the gradient Vx of x with respect 
to X (indicated with F) is such that detF > 0. 

As anticipated above, information on the substructure of each material ele- 
ment are given through the assignment of an order parameter v, by means of 
a sufficiently smooth mapping 89X i — ► v (X) G A4, where A4 is a finite di- 
mensional differentiablc paracompact manifold without boundary endowed with 
metric and connection (Capriz, 1989). The choice of M determines the charac- 
teristic features of each special model of substructure. Vector order parameters 

4 Notations. Some standard notations arc summarized here. Non-standard notations are in- 
troduced in detail in the rest of the paper. Let A and B be tensors of type (p, q) of components 
e.g. A' 1 '"'P . We denote with A • B the standard scalar product given by A' 1 '"'? B* 1 '"'" . In 
particular, if A and B are second order tensors, we denote with AB the product which con- 
tracts only one index and bears a second order tensor; for example, we have (AB) - = A^B^j . 
If A and B arc third order tensors, we indicate with A : B the product contracting two in- 
dices and bearing a second order tensor. If A is a tensor of the type (p,q), with p, q > 0, 
and B is another tensor of the type (r, s), with r, s > and r < p, s < q, or (r = p, s < q) or 
(r < p, s = q), we indicate with AB (with some slight abuse of notation with respect to the 
product between second order tensors) the product which contracts all the indices of B; in 
particular, if p = or q = we take r = and s = respectively. Given two vectors a and 
b, a (g> b denotes their tensor product. In particular, if A and B are second order tensors we 
have AB- (a ® b) = A T a ■ Bb. For any region b of the space, db represents its boundary. 
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with unit length may represent stiff micro-fibers in composites with a softer ma- 
trix or magnetostrictive materials. Vector order parameters not constrained 
to have unit length are used for ferroelectrics, elastic microcracked bodies, ne- 
matic elastomers. Second order tensor valued order parameters may also serve 
as descriptors of families of polymeric chains or polymer stars. 

Let x (X,i) and i> (X,i) be the velocity and the rate of the order parameter 
(in their referential description) evaluated by a given observer. After a change 
in observers ruled by SO (3) new rates x* and i>* are measured: they are given 

by 

x* = x + c(t)+q(t) x (x-x ), (1) 

i>* = v+Aq, (2) 

where c and q are translational and rotational velocities, respectively; .4q is 
the infinitesimal generator of the action of SO (3) on M and A = | q= o (a 
linear operator mapping vectors of R 3 into elements of the tangent space of A4) 
is represented by a matrix with three columns and a number of lines equal to 
dim M, being i/ q the order parameter measured after the action of SO (3). 

2.1 Cracks 

We imagine that the reference place B be free of cracks. When a crack is gen- 
erated in the current configuration of the body, the mapping X i — ► x (X) is 
pointwise one-to-one except a surface C = {X G B, f (X) = 0}, with / a smooth 
function (Fig. 1). The assumption of smoothness for / is only of convenience. 
Notice that C is only a geometrical (non-material) picture in the reference con- 
figuration of the real crack occurring in x (B). 

The intersection of C with the boundary of B is a regular curve OB n C 
endowed with unit normal m such that m (X) belongs to the tangent plane of 
C at each X G dB. Subsets b of B are called 'parts 1 here when they are regular 
regions. When we consider any part be intersecting C, the intersection dbcC\C is 
a regular curve whose normal in the tangent space of C is also indicated with m. 
The normal m to C is defined by m =|^y 5 the opposite of its surface gradient, 
namely L = -Vcm, is the curvature tensor, its trace is the opposite of the 
overall curvature K. In the case treated here, C does not cross completely B. 
The image in B of the real tip of the crack is thus the margin J of C within 
the interior of B. We assume that J is a simple regular curve parametrized by 
arc length s G [0, s] and represented by a point- valued mapping Z : [0,s] — > B 
so that the derivative Z >s (s) of Z with respect to s is the tangent vector t (s) 
at Z(s), while f) = — Z ss is the curvature vector. A normal vector field n is 
chosen along J to be at each Z (s) an element of the tangent plane of C at Z (s) 
outward C (Fig. 1). 

When the crack evolves in the current configuration, its picture in the ref- 
erence configuration is a surface C (t) growing in a certain time interval [0,t]; 
J has then an intrinsic fictitious relative motion with respect to the rest of the 
body, while any piece of C (t) far from J remains at rest. 
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Figure 1: Open crack in the current configuration mapped by inverse motion in 
the reference configuration. 



The intuitive behavior of the crack during the motion is simply described in 
the reference place by the monotonicity of C (t), namely C (t\) C C (t 2 ), Vti < t 2 - 
We assume also that during the time interval in which we study the motion of 
the crack, it does not cut completely the body. In B the velocity of J is 

_ dZ(s,t) 

v tlp - — —. (3) 

Only its normal component V = v t j p • n is independent of the parametrization 
s, and in what follows, we shall consider only v tip = Vn. 

For any field e (X,i) continuous in time and space except C, where it suffers 
bounded discontinuities, its jump [e] there is defined by the difference between 
the outer and the inner trace, i.e., [e] = e + — e~ (when the difference makes 
sense), while the mean value (e) is given by (e) = i (e + + e _ ). When the crack 
is closed in the current configuration, the requirement that its sides do not 
penetrate one into another during the deformation is then [x] • m =0. Some 
special choices of the order parameter require the continuity of it across C. 
Without loss of generality we can consider v continuous across C, while its 
rate may suffer bounded discontinuities. When it is not so, since M is a non- 
linear manifold, the jump u + — v could not make sense and to define [v] it 
should be necessary to embed A4 into an appropriate linear space (a procedure 
based on Withney's or Nash's theorems of embedding). In that case, since the 
embedding is not unique, one should select the one convenient to maintain the 
gauge properties of the underlying physics. 

When there exists any field e (Z,i) defined along J (t) and such that e (X,i) — > 
e (Z (s, t) ,t) as X — > Z (s, t) uniformly in time, we say that e has uniform limit 
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at the tip and confuse e (Z (s, t) ,t) with e (Z (s, t) ,t). In this sense, we indicate 
the tip rate of change of the order parameter as Wti p by considering it at each 
Z Gj as the uniform limit limx^z v (X,i). 

Rates following the crack tip may be defined (by chain rule) as 

= x + Fv hp , i/O = *>+ ( Vi/) v tlp , (4) 

and of course the derivatives x^ and i/O are meant for points away from the 
tip, being rates perceived by observers sitting on the crack tip. Their uniform 
limit at the crack tip are indicated with v tip and w tip respectively, being v tip 
the velocity of the deformed tip and wu p the rate of v> at the deformed tip. 

Moreover, let b be any part of B. The boundary db of b is a two-dimensional 
surface (of normal n) and may be parametrized by parameters u\ and u 2 - If 
we consider b varying in time, i.e., b (t), points X of the boundary 9b (t) are 
identified by X (ui,u 2 , t) so that the velocity u of db (t) is given by 

dX(u 1 ,u 2 ,t) 

u= oi • (5) 

Only the normal component U = u • n is independent of the parametrization 
(til, "2)- Rates following the moving boundary db (t) are then given by 

x°=x + Fu, u° = 0+ (Vip) u. (6) 



3 Balance of standard and substructural inter- 
actions via an invar iance argument 

In the standard mechanics of deformable bodies, common stresses and bulk 
forces are the sole measures of interaction. When the material substructure is 
accounted for, the picture of the interactions become more articulated because 
we must consider objects measuring the extra power developed in the rate i>. 

Below, we summarize balance equations for standard and substructural in- 
teractions in the bulk, at the lateral margins of the crack and at the tip (and 
the attention is focused on the substructural ones because the others are well 
known). To obtain them we could follow two ways: (i) we could use variational 
arguments involving a Lagrangian and its invariance with respect to general 
Lie groups (underlining in this way their covariance 5 ) or (ii) we could require 
SO (3) invariance of the power of 'external' interactions. Though we prefer the 
former, we sketch here the latter because, in this way, we may clearly underline 
the distinction between balance equations and the constitutive structure of the 
interactions involved (two aspects that are mixed in an Hamiltonian approach) . 
Such a distinction clarifies our aim to obtain relations valid even for general 
irreversible processes in deformable solids in different circumstances, although 
our treatment deals mainly with non-linear elastic processes, unless otherwise 
stated. 

5 Following in this way the guidelines of the basic proof of covariance developed for simple 
materials in (Marsden & Hughes, 1983). 
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3.1 Balance of interactions in the bulk 

Let b be any arbitrary part of B (of boundary db) far from the crack. The 
external power V^ xt of the standard and substructural interactions on b is a 
linear real functional on the pairs (x, i>) represented by 



where the substructural interactions are measured through volume (3 and sur- 
face iSn densities as in the case of standard interactions. The bulk density 
[3 may account for both possible substructural inertia effects and interactions 
due, e.g., to electromagnetic fields acting on the substructure (as in the case 
of ferroelectrics) . Both b and (3 are continuous on b. S is called microstress 
and maps linearly vectors of M 3 into elements of the cotangent space of A4. 
Roughly speaking, r = Sn is a 'generalized traction'; the product r ■ v is the 
power exchanged between two adjacent parts at X through a surface of normal 
n, as a consequence of the change of the substructure at the same point. Since 
our analysis is developed in the reference place B, the first Piola-Kirchhoff stress 
T is used: it associates tensions in x (£>) to vectors in B and is the pull-back in 
B of the standard (Cauchy) stress a measuring " true" tensions in the current 
place x(i3) of the body: in fact, T = (detF) <xF~ T . Accordingly, the co-vector 
b is the pull back of the body forces (including inertia) living in the current 
place x (£>). 

We impose that P% xt is invariant under SO (3), for any b, i.e. 



for any choice of c(i), q(t) and b. 

By using (1) and (2), thanks to the arbitrariness of c and q, we obtain the 
standard integral balance of forces 





(8) 




(9) 



and a generalized integral balance of moments: 




(10) 



From (9) the common pointwise balance 



b + DivT = 



in B 



(11) 



follows thanks to the arbitrariness of b, while, from (11), we get 



A T {DivS + f3)= eTF T - (V^ T ) S. 



(12) 
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with e Ricci's alternator. The condition (12) implies that the co- vector qabcTP Ff — 
(V aA%) 6 belongs to the range of the linear operator A T at each v. How- 
ever, A is not one-to-one, then at each v there may exist an element z of the 
cotangent space of M. at u such that 

A T z =eTF T - (V^ T ) S, (13) 

which implies 

DivS - z + f3 = in B. (14) 

Equation (14) is the pointwise balance of substructural interactions (Capriz 
1989) and z is an internal self-force. Equation (13) states that the presence of 
substructural interactions renders unsymmetrical the Cauchy stress <x given by 
er =(dctF) _1 TF T . The case of scalar order parameters seems to be patho- 
logical for the procedure used here because A vanishes. Formally, one may 
circumvent the problem by making use of spherical second order tensors, ob- 
taining the balance (14), then reducing it to the scalar case (which would be in 
any case the Euler-Lagrange equation of some Lagrangian). If one is skeptical 
about such an interpretation, one could accept the point of view discussed here 
to derive (14) only for other types of order parameters and postulate (14) a 
priori in the scalar case. 



3.2 Balance of interactions along the sides of the crack 

Let be be any arbitrary part intersecting the crack C far from the tip. If one 
writes and requires SO (3) invariance, the integral balances (9) and (10) 
follow but be is now the domain of integration. If we shrink be to be H C, since 
the integrands of the volume integrals are continuous while the stress does not, 
we obtain 

f b^O as b c ^ b c n C, (15) 
[ Tin / [T] m as b c -» b c n C. (16) 

Jdbc Jb c nC 

Consequently, the arbitrariness of be implies from (9) the common pointwise 
balance 

[T] m = along C, (17) 
An analogous reasoning can be applied to (10) and leads to 

/ A T [S] n = (18) 

Jdbc 

as be — ► bcHC. The arbitrariness of be implies the pointwise balance A T [S] m = 
which is tantamount to write 

[S] m = z c along C, with A T z' c = 0. (19) 

6 A,B... denote components in B, components in x(B) and a... components over the 

atlas on M. 
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(balance of sub structural interactions along the sides of the crack). z' c vanishes 
when at each u the range of A covers the whole tangent space of M there. 

3.3 Balance of interactions at the tip of the crack 

Take a part of B with the form of a 'curved cylinder' (see, e.g., Fig. 1) b* R 
obtained by translating a disc Dr of diameter R from Z(si) to Z (52) (two 
arbitrary points of J with Si < S2) maintaining Dr orthogonal to t (s) at each 
Z (s) and the centre of Dr coincident with Z (s) at each s. The external power 
of all interactions acting over b* R is 

V^f (A, u, v tlp , w tip ) = f (b • x + f3 • u) + 
+ (Tn • x+Sn ■ v) + I (b tip ■ w Up + Pup-™u P ) , (20) 

J ab* R Jsi 

where h t i P collects only the inertia effects at the tip, and f3 tip accounts formally 
for possible substructural inertia effects at the tip. 

The changes of observers described in (1) and (2) can be written for the 
velocities at the tip as 

v« P = v«„ + c (t) + q (t) x (x 4jp - x ) , (21) 



< p = w tip +A Jp q, (22) 

where „4 t i P = U=o ^ n tne sense °f uniform limit. We impose here SO (3) 

invariance requiring that 

V e b f (±*,v*,v* tip ,w* tip ) = P e b f (x, v, y tip , w« p ) , (23) 

for any choice of c(t), q(i) and b* R . From (23), thanks to the arbitrariness of 
[51,52]) by shrinking b* R to the tip (R — > 0), we obtain 

b tjp + /* Tn = (24) 

tip 

(balance of standard forces at the tip), 

Ptip + / 5n = z« p) with ^ p z; ip = (25) 

J tip 

(balance of substructural interactions at the tip), where we use the notation 
J ti (•) for lim/j^o J D (■)> interpreting it in the sense of uniform limits. 
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3.4 Standard and substructural inertia effects 



3.4.1 Bulk terms 

The bulk interactions b and (3 contain both inertial and non-inertial terms, 
namely b = b OT +b m and (3 = f3 m +(3 nl . The overall power of the inertial terms 
over any stationary part b far from the crack is equal to the opposite of the 
rate of kinetic energy of b, the density of which is commonly chosen to be 

the sum p (j^ ||x|| 2 + k (u, i>)^j , where k(v,i>) the kinetic contribution of the 

substructure, if there is some experimental evidence of it. From the arbitrariness 
of fa, the common identification of the inertial terms follows (Capriz, 1989): 

b in = -px, in = ~pdox + p9 v x, (26) 

where x IS the substructural kinetic co- energy density: its Legendre transform 
with respect to the rate coincides with the kinetic energy k{u,u). In common 
cases, the experiments show that the term jjjpdoX- pdvX 1S m general negligible, 
unless the substructure oscillates at very high frequencies. 



3.4.2 Tip effects 

The standard results before allow us to identify explicitly the tip inertial terms. 
We consider a part around the tip the 'curved cylinder' b* R varying it in time, 
b* R (t), to follow the growth of the crack. We then write the integral versions 
of (26) over b* R adding not only a tip inertial term btip, &s in the common 
treatments, but also a substructural tip inertial term (3 tip , and obtaining then 

/ b" l + r\ Up = -± f p±+ f p iOJ, (27) 
Jb' R (t) J Sl dt h' R (t) Jdb* R (t) 

I ?"+ ["(3 =-± [ pdo X + 
Jb R (t) J Sl at J b , R(t) 

+ I pduX + f pdoxU. (28) 

Jb- R (t) Jdb- R {t) 

The last integral in (27) is the inflow of standard momentum through the bound- 
ary db* R (t) of fajj, due to the 'virtual' (non-material) movement of b* R in time; 
an analogous meaning has the last integral in (28). By shrinking b* R at the tip 
3 (letting R — ► 0) , we get 

b tip = / px (vtip • n) , (29) 

J tip 

Pu P = / pdox »>) (v«p • n) , (30) 

J tip 

thanks to the arbitrariness of [si,S2]- The obvious changes in (24) and (25) 
follow. The result (29) is standard in the theory of dynamic fracture in simple 
bodies. 
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3.5 Stresses at the tip 



The speed of the crack growth is finite, so we may assume that J gD n ® pk be 
bounded up to the tip as R — > 0. Such an assumption (which is in certain sense 
on the behavior of the solution) implies that 



/ 

J til 



n®/9x = 0. (31) 

/ tip 

Moreover, we assume also that the tip flow of substructural momentum be 
bounded up to the tip, i.e., J gD pd^yJJ is bounded as R — > for any choice of 
the order parameter. This implies that d„x lS bounded as R — > (because U is 
bounded) and then 

J n ® pdox = 0. (32) 

J tip 

If the range of Au p covers the whole tangent space of M. at v (Z, t), we get 
z' Up = 0. Consequently, from (24), (25), (29), (30), we get 



Tn = 0, / Sn = 0. (33) 

f tip J tip 

We could reverse the point of view, following a remark of Landau and Lifsits, 
and we could say that since we allow the crack to evolve, the stresses are bounded 
up to the tip. This circumstance together with the assumption (32) would imply 
(33) directly and (31) and z' tip = as further consequences. 



4 Interactions due to the growth of the crack 

The evolution of the crack is represented by the 'fictitious' growth of C in the 
reference configuration B (which on the contrary would remain fixed once and 
for all). Such a growth generates an independent kinematics in B and inter- 
actions power-conjugated with it in the bulk, at the lateral sides and at the 
tip of the crack. They should satisfy appropriate balances. These interactions 
have a twofold nature: from one hand they live in B and are thus different from 
the standard and substructural interactions (notice, e.g., that (11) is the bal- 
ance of forces living in x (£>)); from the other hand, since the kinematics of C 
is only 'apparent' (non-material) in the reference configuration, the new 'forces' 
may be expressed in terms of standard and substructural interactions and the 
free energy. We list below these 'fictitious forces' and their balances that are 
commonly justified in various manners and used to describe different cases of 
mutations in bodies. However, we do not give contribution to the current dis- 
cussion about the attribution of such balances 7 . In the identification of them 
in terms of standard and substructural interactions, our contribution relies in 
the deduction of the substructural components, and this step is crucial toward 
the main result of the present paper. 

7 See the results and the theoretical discussions in (Abcyaratnc & Knowles, 1990; Eshclby, 
1975; Epstein, 2002; Gurtin, 1995; James, 2002; Scgcv, 1996; Silhavy, 1997). 
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4.1 Balance in the bulk of forces due to the crack growth 



For any time varying part b (t) (remind that the evolution of b is not mate- 
rial) we consider bulk internal and external forces, namely the vectors g and 
e respectively, and a stress P, a second order tensor that maps at each X the 
tangent space there onto the corresponding cotangent space. It is commonly 
postulated that they satisfy the integral balance 

/ Pn+ f (g + e) = 0, (34) 

Jdb(t) Jb(t) 

for any choice of b (t) . The pointwise balance 

DivF+g + e = 0, in B, (35) 

then follows. 

In subsequent steps, the strategy foresees the identification of the various 
elements of (35) in terms of the true stresses by making use of a mechanical 
dissipation inequality (a mechanical version of the second law) . We use such a 
procedure having in mind the need to manage an approach valid even in fully 
dissipative processes like viscosity or plastic flows. However, we develop the 
identification of the terms of (35) just in non-linear elasticity (suggesting also 
the necessary developments in other circumstances) to put in evidence the basic 
fact that in that case, the balance (35) reduces to one of the conservation laws 
that can be deduced from a Nother-like theorem, precisely the one associated 
with the invariance of the Lagrangian with respect to diffeomorphisms altering 
B. 

4.2 Balance along the sides of the crack of forces due to 
the crack growth 

For any arbitrary part be (t) intersecting the crack far from the tip, we write 
the balance (34) adding surface measures of interaction along the margins of 
the crack, namely a surface stress u (I — m <g> m) , with I the unit second order 
tensor and a (X) a scalar function continuous up to the tip where it is indicated 
with a u P ] and an internal surface force gc (vector). So that we obtain 

f Pn+ / (g + e)+ f am+ f g c . (36) 

Jdbc(t) Jbc(t) J(dbc(t))nC(t) Jbc(t)nC(t) 

The arbitrariness of be and the bulk balance (35) imply 

[P] m + g c + V c ct + a (O - L) m = along C, (37) 
where Vc denotes the surface gradient (see Section 2). 
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4.3 Balance of tip forces due to the crack growth 



Along the tip we consider a line tension \u p t, a tip internal force g tip (vector), 
and a tip external inertia! force e tip (vector). If we consider a 'curved cylinder' 
b* R (t) (of the type used above) around the tip, intersecting J in two points, 
say Z(si) and Z(s2), in writing on b* R (t) the balance of interactions power- 
conjugated with the 'fictitious' kinematics of C (t), we add a line term 

(Ati P t(s 2 ) - A tip t(si))+ / (gtip+etip)- (38) 

J Sl 

to the bulk and surface terms used in (36). By shrinking b* R (t) up to the tip 
letting R — > (uniformly in time), the arbitrariness of [si,S2] implies 8 

gtip+eup - a Up n - \ tip b, + Pn = along 3. (39) 

J tip 

4. 4 Identification of the inertial term e tip 

To identify explicitly e t i P in terms of the standard and substructural measures 
of interaction, it suffices to consider a 'curved cylinder' b* R (t) wrapped around 
the tip and to write an inertial balance of the type 



rrate 



Jb* R (t) 



rS2 

+ / {b U p ■ v ti p + f3 Up -\v U p + etip ■ vtip) = 0, (40) 

J Sl 

where K ra * e (b* R (t)) is the rate of the kinetic energy of b* R (t) and, of course, 
b t ip and f3 tip develop power in the actual rates. Since b* R varies in time, K ra * e is 

the difference between the time derivative of the integral of p ||x|| 2 + k (u, i>)j 

and its inflow through the moving boundary db* R (t) . By shrinking b* R up to the 
tip (uniformly in time), the arbitrariness of [si,S2] and the use of the identities 
(26), (29), (30) and (32) imply 

e Up = J^p(^\\n 2 + k(v,of)n. (41) 



5 The mechanical dissipation inequality and its 
consequences 

Consider Vu and S. We define the product * by 

(Vv T *S) n • u =Sn • ( Vu) u, (42) 

8 The terms Vcct + ct (KI — L) m and e t i p mark the difference with the equations similar to 
(35), (37) and (39) discussed for the quasi-static evolution of planar three-dimensional cracks 
in (Gurtin & Shvartsman, 1997). 
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for any pair of vectors n and u. Some special cases are the following: when u 
is scalar, \7v T *S = Vu ® S, when (p is a vector, Vv T *S = \7i> T S, while when 
v is a second order tensor, Vv T *S = Vu T : S. 

5.1 The formal statement of the mechanical dissipation 
inequality 

An isothermal version of such an inequality is the following: 

qrate pext < (43) 

where i[<™* e is the rate of the Hclmholtz free energy of b the power of interactions 
over b. The standard assumption is that the (here purely mechanical) inequality 
(43) holds for any choice of the rates involved and for any part b. 

Below we account for time varying parts b (t) of B to follow in Bo the growth 
of the crack in B, then the external power V^ xt must account for the interactions 
associated with such a growth. 

5.2 The mechanical dissipation inequality in the bulk. Con- 
sequences. 

Let b (t) be any time varying part of B far from the crack. is expressed only 
by means of a bulk free energy density ip and we have 

/ V>+ / V(u-n)- / (b-x + /3-i>)- 

Jb(t) Jdb(t) Jb(t) 

- / (Tn • x + Sn ■ i> + (P+F T T+Vv T *S) n • u) < 0. (44) 

Jdb(t) 

where we have used (6) to transform the original integrand Tn • x° + Sn ■ u° + 
Pn • u appearing in the last integral. Since only the component of u normal 
to the surface 9b (t) is independent of the parametrization of db (t), a natural 
invariance requirement with respect to such a parametrization implies that the 
vector (P+F T T+Vi' T *<S) n must be purely normal to 9b; then there exists 
a scalar w such that P+F T T+Vis T *S = zul. By substituting ml within the 
previous inequality, we find the term (ip — w) (u • n). However, since (44) is 
assumed to hold for any choice of the velocity fields, we get tp — w, then 
P = ipl — F T T—\7v T ±S. Notice that, in absence of prominent effects of the 
material substructure described by v, the second order tensor P reduces to the 
well-known Eshelby tensor -01 — F T T. 

We restrict our analysis to the non-homogeneous purely non-linear elas- 
tic case and assume constitutive expressions of the form T = T (X, F, v,Vv) 
for the Piola-Kirchhoff stress, z = z (X, F, for the self-force and S = 

S (X, F, i/,Vi/) for the microstress. If we take the free energy as ip = ip (X, F, i>,Vv), 
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and assume that it admits partial derivatives with respect to its entries, with 
the use of previous results, the mechanical dissipation inequality reduces to 

/ ((0fV - T) ■ F+ (cW - z) • u) + [ {(pdv v il> - 5) ■ Vu) < 0, (45) 
Jb(t) v ' Jb{t) 

where d y ip means partial derivative of ip with respect to the argument y. Its 
validity for any choice of the rates implies 

T = &Fip ; z = ^V ; S = d w „4>. (46) 

As a consequence, taking into account the explicit expression of P, from (35) we 
get g = -obcV and e = -F T b- {Vu T ) (3. 

Special expressions of tft are of Ginzburg-Landau type. In most cases, in fact, 
it seems to be natural to assume ip = ip (X, F, u) + |a (X) || Vf || 2 . In particular, 
if we have V>(X, F,v) = ip 1 (X, F) + xj) 2 (v), » constant and \j) 2 i v ) a coarse- 
grained (perhaps multiwell) energy, the balance of substructural interactions 
(14) coincides with the well known Ginzburg-Landau equation. 

Viscous effects may occur at the gross scale and at the substructural level. 
In this case, the measures of interaction, T, z, S may depend on the rates 
of the fields and their gradients. We assume, as a prototype example, that 
only the self- force z depends on the sole rate i>. The self- force z may be thus 
decomposed into its viscous (v) and non- viscous (nv) parts, namely z = z v +z nv , 
with z nv = z nv (X, F, u,Wu) and z" = z v (X, F, u,Vu; v), with z v i> >0 for any 
choice of i>, which implies z v = A (X, F, is, with A some positive scalar 

function and z nv satisfying (46b). 



5.3 The mechanical dissipation inequality along the sides 
of the crack. Consequences. 

For a part be (t) crossing the crack away the tip, we consider an additional 
surface free energy density 4> (X) along the margins of the crack; it is continuous 
up to the tip where it is indicated with cf> tip . As a consequence, to the bulk terms, 
the ones in (44), we must add the piece 

/ j>+ [ <p (m • u) - / a (m • u) . (47) 

Jb c (t)nc(t) J(dbc(t))nc(t) J(db c (t))nc(t) 

We may then reduce the resulting inequality by shrinking be (t) to C (t) , and 
taking the limit uniformly in time. In this case we get 

/ (<f> - a) (m • u) - [ ([Tm-x] + [5m •!>])< 0. (48) 

The validity of such an inequality for any choice of velocity fields implies (j) = a 
and the local dissipation inequality [Tm • x] + [5m • ii\ < 0. 
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5.4 The mechanical dissipation inequality at the tip of the 
crack. Consequences. 

Let b* R (t) a 'curved cylinder' wrapped around the tip as used in previous sec- 
tions; its boundary intersects the tip in two points Z (s\ (i)) and Z (s 2 (t)). I n 
writing the mechanical dissipation inequality on b* R (t), we consider, in addition 
to bulk and surface energies, a line energy density (; moreover, since the tip 
moves, we must account also for the power of X t i P t and e tip (g t i P is excluded 
because it is internal). Consequently, to bulk and surface contributions we add 
the term 

*2(t) \ 

C ~ (Atipt • v t i p \ S=S2 — Xti P t ■ v t ip | s = Sl ) . (49) 
«i(*) / 

By shrinking b* R (t) up to the tip, taking the limit uniformly in time and making 
use of the line balance (39), we get 

(C — X tip ) (t • V t i p ls=s 2 ~ * ' Vfip |s=si ) + 

+ / Cf)-v tl p + / v Up ■ g Up < (50) 

Since the resulting inequality must be valid for any choice of the velocity fields, 
then of Vtip, we obtain £ = X t i P and gtip'^Up < 0, which reduces to Vgu P • n < 
since vu p = Vn, implying then that the component of g t j p along the motion of 
J, namely g tip = g t i P • n, must have a structure of the type gu P — Qu P V, where 
Qu p is a negative 'diffusion' coefficient that must be assigned constitutively. 

6 Driving the tip of the crack 
6.1 The driving force 

The explicit expressions of P, g, a, Xu p allow us to write the tip balance (39) 



as 



-<j> Up n - X tip t)+ jf (pQ||i|| 2 + /.-(V.i>))l-P]n = - g„„. (r,l) 



We indicate with j the vector 



Lip 



p[h\±\\ 2 + k(v,u))I-F)n (52) 



It represents the tip traction exerted by the bulk material on an infinitesimal 
neighborhood around the tip. Let n (s) be the direction of propagation of the 
crack at the point Z (s) of the tip, the component of (51) along n (s) is given by 

-Kp ~ x ti P & + n • j = -QtipV. (53) 
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By indicating with J (the J-integral) the product n j, we interpret the difference 

f = J - 4>u P - (54) 

as the force driving the tip of the crack; it accounts directly for the influence of 
the material substructure. Since g tip • v t ip < 0. we get fV > (wc remind that 
•vup is of the form Vn). When the crack grows, i.e., when V > 0, the driving 
force must be non-negative, i.e., f > 0. 

6.2 Dynamic energy release rate at the tip. 

The energy release rate at the tip is given by the power fV developed by the 
driving force along the normal motion of the crack tip. Previous results allow us 
to express the product fV in terms of the power of standard and substructural 
interactions and of the free energy, and we obtain 

-4> Up V - \ tip SW = ft/, (55) 
which represents the balance of energy at the tip. 

6.3 Quasi-static extended J-integral and its path indepen- 
dence 

When inertial effects are negligible, the evolution of the crack is "quasi static" . 
The J-integral (J = n • j) reduces to its quasi static counterpart J 9 „: 

J gs = n • f Pn. (56) 

J tip 

It reduces to the standard J-integral given by n • J Up (tJA — F T Tj n when 

the substructure is absent or its gross effects are negligible. 

If the material is homogeneous, C is planar (i.e. the crack is straight), the 
crack has the margins free of standard and substructural tractions ( in the sense 
that T ± m = and S ± m = 0), J qs is path-independent. 

The hypotheses of homogeneity of the material and absence of inertial effects 
imply g = and e = 0. Then, the bulk balance (35) reduces to DivP = 0. 
Moreover, since the margins of the crack are free of standard and substructural 
tractions, we get [P] m = [ip] m. 

With these premises, we take in B an arbitrary 'curved cylinder' b* R wrapped 
around the tip where we have 
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Dr is arbitrary, then we need only to evaluate the difference 

n (s) • j Pn-n (s) • f Pn, (58) 

J dD r J tip 

that is equal to n- J DRnC ["0] m ) which vanishes because the crack is straight, i.e. 
m ±n. The path-independence of ^ qs follows. 

7 The energy dissipated into a process zone of 
finite size around the crack tip 

When a crack propagates, a material part Pz around the tip becomes highly 
unstable, in certain sense 'fragmented' (Aoki, Kishimoto & Sakata, 1981; 1984). 
Usually, Pz is called process zone. In brittle fracture, the process zone may be 
considered practically coincident with the tip, while in ductile fracture Pz has 
finite size. For the latter case we obtain new path-integrals which allow us to 
evaluate the energy dissipated during the evolution of the crack. 

We assume that Pz can be approximated reasonably by a 'curved cylinder' 
P wrapped around J (basically, P has geometrical properties analogous to b* R 
used in previous sections); the intersection of P with the plane orthogonal to 
the tangent t (s) of J at s G [0, s] is a disc P*' with the centre on J (being the 
centre the sole intersection of J with P). The approximation of Pz with P is 
rather rough but it does not influence the basic structure of the results obtained 
in the present section. During the evolution of the crack, P varies in time and is 
P(t). The boundary dP is a surface (with outward unit normal indicated with 
n) parametrized by v\, vi and points ~X.(vi,V2,t) G dP (t) have an intrinsic 
velocity u given by u = <9tX (i>i, V2, t). Consequently, rates following dP (t) are 
x > = x + Fu, and = j>+ (Vis) u. 

The velocity u at the boundary dP is decomposed as 

u = u tr (t) + q (t) x (X - X ) + a (t) (X - X ) + u d (X, t) , (59) 

Utr (t) denotes the component of rigid translation; q(t) x (X — X ) the rota- 
tional component, with Xo an arbitrary fixed point (note that the presence 
of (X — Xo) instead of (x — xo) is due to the circumstance that u is a mate- 
rial velocity in B); a (t) (X — X ) the velocity associated with the self-similar 
expansion of P; (X, t) the component of the velocity due to the distortion. 

Let now bjj be another 'curved cylinder' wrapped around the tip, fixed in 
time and containing a piece P of the process zone. We analyze the behavior of 
P (t) in a time interval in which dP (t) does not intersect db * R . 

With reference to the situation described before, we indicate with <3E> (P) the 
rate of energy dissipated in P during the evolution of the crack and assume that 
all the dissipation is concentrated in P during the evolution of the crack. Such 
an assumption implies that 

i>= [ ( T • F + z • u+S ■ ViA . (60) 
* R \P(t) Jb R \P(t) v > 
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In other words, the mechanical dissipation inequality in b* R \P (t) reduces to an 
equation because no dissipation mechanism occur outside P. 

With these premises, the balance of the energy takes the form 



J t ^ m P Q \\n 2 + k (v, u)j + f g _^ p Q ||x|| 2 + k (i/, *)) (u • n) + 

+^ / i>+ ! V(u-n)-$(P) = 

at Jb* R \P(t) JdP(t) 

= f (b" 4 -x + (3 m -0) + f (Tn • x+5n • u) . (61) 

Jb R \P(t) Jdb R 

By making use of the weak form over b^\P (t) of the balances of standard 
and substructural interactions (11) and (14), and of the integral identity 

/ (Tn • x+5n ■ if) = / (Tn • + Sn ■ i/>) + 

Jdb* R JdP(t) 

+ [ {Div(±T) + Div{vS)), (62) 
Jb R \P(t) 

from the arbitrariness of the piece of J considered, we obtain 



where 



$ (P T ) = u tr (t) • j (P T ) + q (i) L + a (t) M + I, (63) 



l(P")= I [p[l\\n 2 + k(u,0)\l-P\n ((,4) 



dP"{t) 



1 = ^ t (X - X ) x (p Q ||xf + fc (i/, I-p) n (65) 
M =J dp t (p(l\\xf + k(v,v)^)I-^)n-(X-Xo) (66) 



1= 



IdPT(t) 



n • u d 



dPT(t) 



(Tn • x + <Sn ■ if) . 



(67) 



8 Special cases 

The theory discussed in previous sections allows us to describe the behavior of 
cracks in several cases of complex materials. Two essential ingredients are nec- 
essary to apply the results: (i) the choice of an order parameter i> (hence of M) 
describing appropriately the material substructure; (ii) an explicit expression of 
the free energy. In what follows, we indicate two possible spheres of application: 
ferroelectric solids and solids exhibiting strain gradient effects. 
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8.1 Cracks in ferroelectrics 



8.1.1 Preliminary remarks 

The predictions of Griffith's theory fall for cracks propagating in ferroelectric 
solids subjected to the action of external electric fields: there is a discrepancy 
between the driving force predicted theoretically and the experimental data. 
Basically, we regard this discrepancy as due to the circumstance that Griffith's 
theory does not account for the substructural interactions associated with the 
occurrence of spontaneous polarization that can be induced by strain, variation 
of temperature and applied electric fields. The polarization is indicated with p 
and we take as order parameter the vector p = p _1 p such that < |p| < p m , 
with p m a material constant. Then M is the ball of radius p m in M 3 . We also 
consider the body as subjected to an external electric field (S. 

Balance equations are formally identical to (11), (13) and (14). In this case, 
A is the second order tensor px (in components (px)-^ = e^pfe, with the 
alternating symbol). As a consequence, the relation (5) becomes p* = p + pxq. 
Here, the microstress S accounts for interactions between neighboring crystals 
with different polarizations; z measures self-interactions within each polarized 
crystal. To include the effects of the applied electric field in the balance equa- 
tions we assume that the bulk interactions b and (3 and the boundary 'tractions' 
t = Tn and r =Sn can be decomposed additively in electromechanical parts 
(em) and purely electric parts (el): 

b = b em +b e/ ; (3 = f3 em +(3 el ; (68) 



t — t em +t e i ; T T em +T e ; . (69) 

For any arbitrary part b of B, the purely electric parts are characterized by 
the balance 

^-D(b)+ /(b er x + /3 er p)+ / (t er x + r erP ) = 0, (70) 

at Jb Jdb 

where D (b) = \ J b \€\ = —\ J b p£-p is the electric energy of b. 

A theorem of Tiersten (1964) give us the explicit expression of the rate of 
D in the current configuration. We indicate the material version of Tiersten's 
formula by pulling it back in the reference configuration: 

j f D (b) = - J p (grade) p.±-J± (det F)p^F~ T n -i-JpiE-p, (71) 

where p n is the normal component of p, grad indicates the gradient with respect 
to x. It follows that 

b el =p (grad<E) p ; f3 el = p<5; (72) 
t ei = i(detF)p2 F - T n ; r el = 0. (73) 
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Consequently, the balance equations (11) and (14) become 
b em + DivT + p (grad€) p = 0, 



(74) 



(3 em - z+DivS + P £ = 0, (75) 
where b em and (3 em include the inertial terms as ever (see also Davi, 2001). 

8.1.2 The driving force in ferroelectrics 

Taking in mind the identification of u with p, all the general results on cracks 
presented before apply and for non-linear elastic ferroelectrics the quasi-static 
J-integral is given by 

J 9S - n • f Ul - F T G) F ^-Vp T a V p^) n (76) 
Jtip v ' 

An expression of this type has been used in (Wei & Hutchinson, 1997) for the 
special case of infinitesimal strains and fits reasonably the experimental data 
about the driving force. 

8.2 Cracks in materials with strain gradient effects 

Size effects are well recognized in the behavior of crystalline solids even during 
phase transitions. To describe these experimental evidences, models involving 
the second gradient of deformation have been considered (Dunn & Scrrin, 1985). 
Strain gradient effects are induced by a latent substructure. Following the ther- 
modynamically consistent theory of Capriz (1985), latence is induced here by an 
internal constraint obtained by identifying v with F, in absence of external ac- 
tions on the substructure (/3 = 0). As a consequence, if one assumes constitutive 
equations of the type T = f (F, VF), S = S (F, VF) and ip = ip (F, VF), it fol- 
lows that T = <9fV> $ = c'vfV' and the balances of standard and substructural 
interactions merge one into the other and reduce to the sole balance 

Div {T-DivS) + b = 0. (77) 

The quasi static J-integral J gs takes the special form 

J gs = n ■ f Ul - F T a F ^-VF T : n; (78) 

Jtip v ' 

its reduced version in infinitesimal deformations has been used in (Fulton & 
Gao, 2001) to calculate the driving force in crystalline materials that fail by 
decohesion at atomic scale: the results fit reasonably experimental data. 
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9 Discussion 



We have presented a fully non-linear three-dimensional description of the crack 
growth under large strains in complex bodies suffering a prominent influence of 
the material substructure on the gross behavior. We have used 'abstract' order 
parameter fields as coarse-grained descriptors of the morphology of the sub- 
structure. From one hand such a point of view allows us to unify some existing 
preliminary tentative to modify the standard theory of fracture in special cases 
of complex bodies, while, on the other hand, it furnishes a general tool able to 
provide directly the expression of the driving force in all cases of materials ad- 
mitting Ginzburg-Landau-like energies. The present paper extends and renders 
more perspicuous preliminary two-dimensional results presented in a section of 
the article (Mariano, 2001). 
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